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Abstract 


A  minimum  correction  homotopy  approach  is  used  to  obtain  the 
simultaneous/integrated  optimal  design  of  a  large  flexible  structure  and  its  active  control 
system.  Instead  of  the  usual  method  of  weighting  and  summing  all  desired  objectives  to 
form  a  constrained  scalar  optimization  problem,  a  vector  of  objective  functions  is  dealt  with 
directly. 

The  Draper/RPL  configuration  (a  central  hub  with  four  symmetric,  identical  arms)  is 
the  design  structure.  The  design  seeks  to  minimize  the  mass  of  the  arms.  Using  simple 
feedback  of  arm  displacements  and  velocities,  the  control  system  seeks  to  achieve  specified 
closed-loop  eigenvalues  (frequencies  and  damping  ratios)  and  control  effort.  Design 
variables  are  the  arm  dimensions,  control  system  gains,  and  sensor  and  actuator  locations. 
Not  only  can  the  structural  design  be  accomplished  while  placing  the  closed-loop 
eigenvalues,  but  a  simultaneous  50%  reduction  in  mass  and/or  control  effort  can  be 
obtained. 

Since  reduced-order  models  were  used  for  the  structural/control  design,  the  resultant 
configurations  are  easily  driven  unstable  by  spillover  from  higher-order  unmodeled  modes. 
A  modal  suppression  technique  is  applied  to  eliminate  observation  spillover  and  provide  a 
decade  of  deadband  above  the  controller  bandwidth. 
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INTEGRATED  STRUCTURAL/CONTROL  DESIGN 
VIA  MULTIOBJECnVE  OPTIMIZATION 


I.  Introduction 


The  design  of  large  flexible  space  structures  presents  distinct  and  often  competing 
challenges  to  both  structural  dynamicists  and  control  engineers.  The  need  to  maneuver 
such  structures  and  then  suppress  resultant  vibrations,  or  to  reject  disturbance  vibrations 
from  a  variety  of  sources  to  meet  exacting  pointing  and  stability  requirements  is  often 
directly  opposed  to  the  low  stiffness  and  high  flexibility  accompanying  a  low  weight 
design.  Passive  techniques  to  supplement  structural  damping  justly  receive  much  attention 
and  are  very  important,  but  are  not  panaceas.  The  requirement  for  active  control  of  large 
space  structures  is  widely,  if  not  universally,  accepted. 

Active  control  of  large  space  structures  has  been  extensively  studied.  Design 
methods  have  progressed  well  beyond  the  practice  of  adding  active  controls  to  make  up  for 
or  'fix'  structural  difficulties.  However,  design  of  actively  controlled  structures  has 
traditionally  been  a  sequential  process:  first  the  structure  was  designed  based  on  structural 
criteria  (minimize  mass  while  achieving  desired  natural  frequencies,  mode  shapes  and 
dynamic  response);  then  the  control  system  was  designed  to  meet  desired  closed-loop 
objectives  (while  minimizing  controller  work  and  the  energy  of  the  vibrating  structure)  for 
the  given  structural  design.  Inasmuch  as  a  lower  weight  (higher  flexibility)  structure 
requires  more  control  energy,  the  objectives  of  the  two  design  steps  are  contradictory  and 
an  optimal  controller  placed  on  an  optimally  designed  structure  does  not  result  in  an  optimal 
control-structure  design. 


A  natural  progression  has  been  to  incorporate  the  sequential  approach  into  a  design 
loop,  repeating  the  process  so  that  the  design  result  improves  with  each  iteration. 
However,  the  iterated  design  tends  to  be  dominated  by  the  first  discipline  in  the  sequence. 

In  the  recent  technical  literature,  the  competitive  natures  of  structural  and  active 
control  system  design  are  recognized.  Attempts  are  made  and  interfaces  are  suggested  to 
integrate  the  design  process.  However,  because  of  the  very  strong,  frequently 
unintentional  and  adverse  coupling  between  the  flexible  structure  and  active  control  system, 
a  wholistic  or  simultaneous  approach  to  integrated  structural/control  design  is  necessary. 
In  particular,  a  vector  approach  wherein  both  structural  and  control  design  objectives  are 
treated  in  parallel  will  be  developed. 

Finally,  it  must  be  noted  that  designs  are  normally  based  on  reduced-order  models, 
especially  of  large  flexible  space  structures,  either  due  to  considerations  of  on-board 
computer  speed  or  the  ability  to  accurately  determine  controller  gains.  Spillover  from  the 
higher-order  (unmodeled)  residual  modes  affects  system  performance  and  may  well  cause 
instability.  Considerations  of  model  order  truncation  effects  must  not  be  neglected  in  any 
structural/control  design.  An  integrated  design  approach,  as  any  other,  is  only  as  good  as 
the  information  it  is  based  upon. 

This  dissertation  presents  an  approach  to  integrate  structural  and  control  design. 
First,  previous  work,  using  both  scalar  and  multiobjective  approaches,  is  reviewed 
(Chapter  II).  Next,  the  problem  of  an  actively  controlled  structural  system  is  formulated 
(Chapter  III),  followed  by  descriptions  of  how  the  structure  and  controller  are  modeled  for 
this  study  (Chapter  IV).  In  Chapter  V  the  use  of  a  minimum  correction  homotopy 
algorithm  to  optimize  a  vector  objective  function  is  discussed  and  applied  to  this  integrated 
structural/control  design  problem.  Chapter  VI  presents  design  results  for  various 
combinations  of  design  objectives.  An  examination  of  the  designs’  sensitivity  to  spillover 
from  higher-order  modes  leads  to  the  two  methods  of  spillover  control  and  corresponding 
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results  discussed  in  Chapter  VII.  Conclusions  regarding  this  integrated  design  approach, 
including  spillover  control,  and  recommendations  for  future  work  are  given  in  Chapters 
Vm  and  IX,  respectively. 


Within  the  last  several  years,  interest  and  research  in  integrated  structural/control 
design  has  grown  and  been  reported.  Such  design  methods  result  in  desired  system 
response  characteristics  with  lower  mass  and  possibly  control  effort  because  structural 
properties  have  been  tailored  to  augment  the  active  control  system.  However,  the  vast 
majority  of  the  work  has  used  a  series/sequential  approach  to  arrive  at  an  optimal  design. 


Scalar  Approach 

This  integrated  design  task  has  almost  exclusively  been  treated  as  a  scalar 
optimization  problem.  The  objective  function  is  usually  taken  to  be  either:  1)  the  structural 
mass,  with  the  other  desired  criteria  (controller  performance  or  system  response 
characteristics)  treated  as  constraints  to  define  a  feasible  design  space,  or  2)  a  weighted 
sum  of  all  desired  properties  and  characteristics.  In  the  first  case,  such  an  approach  not 
only  requires  a  priori  determination  of  the  criteria  excluded  from  the  objective  function,  but 
seriously  weakens  the  overall  design  process  by  the  choice  of  a  single  criterion  to  define 
the  merit  of  an  entire  system.  Design  trade-offs  may  be  accomplished  (with  some 
difficulty)  by  changing  the  constraints  to  define  a  new  feasible  design  space  and 
reaccomplishing  the  optimization  problem. 

In  the  second  case,  the  combination  of  several  conflicting,  usually  non- 
commensurable  criteria  which  should  be  optimized  simultaneously  into  a  single  scalar 
objective  function  is  not  only  unnatural  to  the  physics  of  the  problem  but  also  often 
inadequate.  While  commendable  and  significant  in  that  structural  and  control  syntheses  are 
treated  together  vice  sequentially,  this  combinatorial  approach  raises  the  question:  what 
does  minimization  of  a  sum  of  structural  properties  (say  mass)  and  controller  characteristics 


(say  control  energy  or  a  steady-state  linear  quadratic  regulator  [LQR]  cost  function)  really 
mean  or  represent?  LQR  controller  synthesis  suffers  a  similar  identity  crisis:  why  the 
minimization  of  a  particular  weighted  integral  of  state  and  control  variables?  In  the  case  of 
LQR  controllers,  the  question  is  moot  since  the  true  motivation  is  to  yield  a  tractable 
problem  whose  solution  is  readily  synthesized  and  easily  implemented.  Likewise,  the 
combinatorial  approaches  to  structural/control  design  "represent  at  least  a  convenient 
parameterization  of  the  problem  wherein  designs  can  be  iteratively  considered  and 
improved  through  variation  of  the  weighting  matrices  (1:1 124)."  Although  doing  such 
design  sensitivity  trade-offs  only  yields  local  information  in  the  neighborhood  of  the 
optimum  (2:483,503;  3:141;  4:1101;  5:333),  this  approach  is  particularly  and  conveniently 
amenable  to  solution  via  the  vast  software  libraries  of  existing  optimization  techniques  and 
algorithms. 


Multiobiective  Optimization 

Multiobjective  (multicriteria,  multicriterion,  vector)  optimization  is  a  more  natural 
and,  hopefully,  more  efficient  approach  to  effectively  account  for  the  numerous  different 
and  often  conflicting  or  competing  criteria  inherent  in  structural/control  design  problems. 
This  approach  may  be  mathematically  considered  as  the  vector  extension  of  scalar 
optimization  or,  from  an  engineering  point  of  view,  seen  as  a  tool  to  find  compromise 
designs/solutions  to  the  conflicting  practical  requirements.  The  designer  can  then 
systematically  analyze  the  alternatives  to  arrive  at  a  preferred  solution  which,  while  none  of 
the  criteria  will  necessarily  attain  its  extremum,  will  satisfy  the  design  requirements  in  some 
subjective  'best'  way  (2:503;  3:141;  4:1101). 

Multiobjective  optimization  first  arose  in  the  study  of  mathematical  economics  and 
progressed  to  general  decision-making  and  engineering  (6:162-163;  7;  8;  9;  10;  11:1-10). 
Its  use  in  structural  mechanics  dates  back  to  the  1970's.  Stadler's  survey  (12)  examines 
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multicriteria  optimization  in  the  broader  field  of  mechanics  to  1984.  Briefer  summaries  of 
its  applications  in  structural  mechanics  to  date  are  given  by  several  authors  (2:484;  3:141- 
142;  13:184;  14:119-120). 


Pareto  Optimality.  (2:484-487;  3:142;  4:1102-1103;  5:333-334;  10:68-72;  15:925- 
926;  16:459-460).  The  general  form  of  a  multiobjective  optimization  problem  is 

min  (2-1) 

ze  Q 

where  y:  Q  — >  9lP  is  a  vector  objective  function  given  by 

:y  (z)  =  [  Yi  (z),  Ya  (z),  •  •  • ,  Yp  (z),  f  (2-2) 

The  components  Yi-  >  51,  i  =  1,  2,  ....  p  are  the  design  criteria  -  the  conflicting  and 
often  non-cotnmensurable  performance  objectives.  The  design  variable  vector  g  belongs  to 
the  feasible  set  Q  <=  defined  by 

a  =  {  z  e  91^^ :  £  (z)  $  0,  li  (z)  =  0  }  (2-3) 

where  Cj:  91^  — »  91^  i  =  1,2, ...,  r  are  inequality  constraint  functions  and  hj:  914  i  - 
1,  2, ...,  q  are  equality  constraint  functions.  TTie  image  of  the  feasible  set  in  the  criterion 
space  is  the  attainable  set,  given  by 

©  =  {Y(z)e91*’:ze  Q}  (2-4) 

Since  the  components  Yi  of  the  objective  vector  are  usually  conflicting,  an  optimal 
solution  (or  superior  solution),  that  is,  a  unique  point  or  value  of  z  which  will  minimize  all 
Yi  simultaneously,  will  not,  in  general,  exist  for  the  multiobjective  problem.  Attention  is 
therefore  commonly  directed  to  Pareto  optimum  solutions  (non-inferior  solutions,  non- 
dominated  solutions,  efficient  solutions). 


A  vector  z*  €  is  (strongly)  Pareto  optimal  if  and  only  if  there  exists  no  £  e  Q  such 
that  Yj  (^  <  Yj  (z*),  i  =  1,  2, ....  p  with  Yj  (z)  <  Yj  (x*)  for  at  least  one  j.  Verbally,  z*  is 
Pareto  optimal  if  there  exists  no  feasible  vector  z  which  would  decrease  some  criterion 
without  causing  a  simultaneous  increase  in  at  least  one  other  criterion.  In  contrast,  if  all 
elements  of  the  objective  vector  are  simultaneously  achievable,  the  solution  is  termed 
"dominated”. 

If  0  is  non-empty  closed  and  max  (Yi  (z) :  Yi  (l)  ^  0)  <  <»  for  all  i,  then  Q  has  at 
least  one  (strongly)  Pareto  optimal  solution.  "Thus,  a  large  class  of  multiobjective 
optimization  problems  in  structural  design  may  be  expected  to  possess  at  least  one  non- 
dominated  solution  (16:460)."  Whereas  in  scalar  optimization  one  optimal  solution  is 
usually  characteristic  of  a  problem,  a  set  or  family  of  Pareto  optimal  solutions  generally 
exist  for  a  multiobjective  optimization  problem.  Mathematically  speaking,  the 
multiobjective  optimization  problem  is  'solved'  once  the  Pareto  optimal  set  is  determined. 
Practically  speaking,  too  many  Pareto  optimal  solutions  often  exist,  and  it  may  be 
necessary  to  order  or  rank  the  set  to  determine  a  preferred  solution. 

While  there  is  an  extensive  literature  of  models  and  methods  to  generate  the  Pareto 
optimal  set  and  hence  solve  multiobjective  optimization  problems,  the  choice  of  which 
technique  to  use  on  a  given  problem  is  very  subjective,  especially  for  nonlinear,  large 
dimensioned  problems.  There  is  no  requirement  to  seek  Pareto  over  dominated  solutions; 
it  has  simply  become  common  practice  to  do  so  when  possible.  If  they  can  be  found,  the 
Pareto  solutions  can  be  claimed  to  be  optimal  in  some  subjective  sense,  whereas  dominated 
solutions  are  only  attainable.  Unfonunately,  there  are  very  few  direct  precedents  for  the 
multiobjective  optimization  of  large,  actively  controlled  structures.  Based  on  the  work  and 
conclusions  of  Rao  and  associates  (13,  14,  17),  utility  function,  goal  programming  and 
cooperative  game  theory  methods  appeared  promising. 


The  utility  function  and  goal  programming  methods  both  reduce  the  original 
multiobjective  problem  to  a  scalar  optimization  problem  (primarily  as  a  weighted  sum  of  the 
individual  design  objectives).  Minimizing  the  scalar  objective  naturally  yields  a  Pareto 
solution  (corresponding  to  the  particular  scalarization  scheme,  weighting,  etc.).  However, 
both  of  these  methods  had  difficulty  converging  to  a  solution  for  the  sub-problem  of 
control  system  design  (placing  closed-loop  eigenvalues)  and  would  not  converge  to 
solutions  for  integrated  structural/control  design  (reducing  mass  while  placing  closed-loop 
eigenvalues).  Therefore,  these  relatively  simple  approaches  were  abandoned  and,  based  on 
these  experiences,  cooperative  game  theory  (a  much  more  complicated  and  elegant 
scalarization  scheme)  was  never  attempted. 

Minimum  Correction  Homotonv  Approach.  Junkins  and  associates  have  used  a 
minimum  correction  homotopy  (MCH)  approach  to  multiobjective  optimization  which  can 
preserve  the  vector  nature  of  the  competing  design  criteria  (18-21).  However,  in  three  of 
the  four  works  (18-20),  only  controller  design  objectives  (state  error  energy,  control 
energy,  and  a  stability  robustness  measure)  were  considered  in  a  sequential  approach 
minimizing  one  criterion  while  holding  the  other  two  objectives  to  prescribed  values  in  the 
neighborhood  of  their  unconstrained  minima.  In  the  fourth  work  (21),  the  MCH  algorithm 
was  used  to  establish  a  feasible  design  point  satisfying  a  constraint  vector  of  desired 
closed-loop  eigenvalues.  A  robustness  measure  (the  sensitivity  of  the  closed-loop 
eigenvalues  with  respect  to  variation  of  uncertain  system  parameters)  was  subsequently 
minimized. 

While  the  most  notable  and  honestly  vector  approaches  to  multiobjective 
optimization,  these  four  studies  (18-21)  still  lack  two  important  considerations  from  the 
standpoint  of  integrated  structural/control  design.  First,  only  control  system  design  was 
considered,  although  one  case  is  presented  (21)  where  a  structural  design  iteration  was 
performed  before  the  MCH  approach  was  used  to  satisfy  the  constraint  vector.  Second, 


the  'true'  objectives  are  not  dealt  with  directly  as  a  vector.  They  are  either  dealt  with  in 
sequence  as  constraints  (18-20)  or  'after  the  fact'  (21).  Indeed,  values  for  the  final  or 
'true'  objective  (robustness)  are  not  even  reponed  in  (21). 

Given  the  convergence  problems  encountered  with  utility  function  and  goal 
programming  approaches  to  finding  Pareto  optimal  solutions,  the  MCH  approach  was 
turned  to  next.  A  simple  extension  of  the  vector  (objective)  function  in  (21)  to  include 
structural  mass  (and  later  control  effort)  applied  the  MCH  technique  to  the  problem  of 
integrated  structural/control  design. 

Although  Stadler  reported  "the  most  extensive  use  of  multicriteria  optimization  has 
been  made  in  optimal  structural  design  (12:282),"  the  techniques  encompassed  by  the  topic 
of  multiobjective  optimization  are  only  beginning  to  be  considered  for  the  simultaneous  or 
integrated  structural/control  design  problem.  While  again  commendable  unifications  of 
structural  and  control  design  objectives,  the  work  to  date  may  still  be  classified  as  either  a 
series/sequential  approach  (17-20)  or  reduction  (via  some  weighting  scheme)  to  a  scalar 
optimization  problem  (13,  14, 17-20). 

The  MCH  method  differs  from  the  traditional  or  standard  multiobjective  optimization 
techniques  in  two  ways.  First,  dominated  vice  Pareto  solutions  are  found;  there  is  no 
trade-off  amongst  the  competing  design  objectives  to  arrive  at  some  "best"  solution.  The 
method  converges  rapidly  to  the  objective  vector  as  specified.  Second,  the  vector  nature  of 
the  objectives  is  preserved  vice  the  more  common  reduction  to  a  scalar  optimization 
problem. 

In  spite  of  these  differences,  the  topic  of  multiobjective  optimization  should  include 
the  MCH  approach  in  its  catalog  of  techniques.  This  research  shows  the  approach's 
validity  and  great  utility  to  extend  multiobjective  optimization  to  the  problem  of  integrated 
structural/control  design. 


III.  Problem  Formulation 


Consider  the  discretized  linearized  equations  of  vibrational  motion  for  a  controlled 


structural  system,  neglecting  disturbances: 


(3-1) 


where 


Mw  +  Cw  +  Kw  =  Bu 


n- vector  of  generalized  coordinates 

a-vector  of  control  inputs 

nxn  symmetric  positive  definite  mass  matrix 

nxn  symmetric  positive  semidefinite  damping  matrix 

nxn  symmetric  positive  semidefinite  stiffness  matrix 

nxa  control  influence  matrix 

d(  )/dt 


For  direct  output  feedback  control  of  such  a  system,  let  local  position  and  velocity 


measurements,  respectively,  be  denoted  as: 


Jlp  =  Hp  w  ,  jiv  =  Hv  w 


(3-2) 


Assuming  s  colocated  position  and  velocity  sensors,  Hp  =  Hy  =  H  and  is  an  sxn  matrix 
describing  sensor  locations  and  orientations,  while  yp  and  y^  are  s- vectors. 

The  control  u  is  a  linear  combination  of  the  outputs: 

u  =  —  I  Gp  yp  +  Gy  yy  ] 

=  -GpHw-GyHw 


(3-3) 


where  Gp  and  Gy  are  constant  axs  gain  matrices. 
If  proportional  damping  is  assumed. 


C  =  a  M  +  p  K  (3-4) 

where  a  and  P  are  positive  constants.  If  Eqs  (3-3)  and  (3-4)  are  substituted  into  Eq  (3-1), 
the  closed-loop  system  is 

Mciw  -I-  Cciw  -t-  K^jw  =  0  (3-5) 


where 

Mci  =  M  ,  Cci  =  a  M  +  p  K  +  B  Gy  H  ,  =  K  +  B  Gp  H  (3-6) 

The  second-order  system  of  Eq  (3-5)  can  be  written  in  first-order  state-space  form: 

Ax  =  Tx  (3-7) 

where 
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or,  more  conveniently,  as 


X  =  A(.|X 


(3-9) 


where,  since  M^i  =  M  and  with  I  an  identity  matrix, 

0  I 

I 

AqI  — 


Ke, 


-M  ’  C, 


cl 


(3-10) 


The  associated  right  and  left  eigenvalues  and  eigenvectors  of  A^-i  are,  respectively, 
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£i  -i  ~ 


i  =  1,2,  2n 


(3-11) 


Assuming  all  eigenvalues  are  distinct  and,  therefore,  the  eigenvectors  are  linearly 
independent,  the  eigenvectors  are  conventionally  normalized  as 

=  Sy,  fi'^Acifij  =  5y>.i  (3- 12a) 

or 

FT  E  =  I ,  PT  Aci  E  =  r  A  J  (3-12b) 

where  E  and  F  are  the  right  and  left  eigenvector  matrices  whose  ith  columns  are  the 
eigenvectors  and  fj,  respectively,  6y  is  the  Kronecker  delta,  and  f  A  J  is  a  diagonal 
matrix  of  the  eigenvalues  Xj,  ^2, ...,  ^211  • 

The  closed-loop  damping  ratios  and  damped  frequencies  cOj  are  related  to  the  in 
general  complex  conjugate  pairs  of  closed-loop  eigenvalues  Xj  as 

Ci  =  -  Cj  /  (  +  0)  2  )l/2 ,  =  a,  ±  j  COj  i  =  1,  2,  ...,  n  (3-13) 

The  specific  structural  model  and  control  system  used  in  this  research  arc  described  in 
the  next  chapter,  leading  to  explicit  expressions  for  the  system  matrices  and  state,  output 
and  control  vectors. 
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IV.  Model  Description 


Integrated  structural/control  design  modifications  to  the  Draper/RPL  spacecraft  model 
will  be  demonstrated  herein.  This  model  has  been  used  to  test  and  demonstrate  hardware 
and  control  laws  for  maneuvering  large  flexible  spacecraft.  As  shown  in  Figure  4-1,  the 
configuration  consists  of  a  large  central  hub  with  four  identical  appendages/arms 
symmetrically  cantilevered  from  the  hub.  Each  arm  is  modeled  as  a  continuous  beam  with 
a  lumped  mass  at  the  end. 


Figure  4-1.  The  Draper/RPL  configuration 


Only  the  planar  rotational  and  vibrational  dynamics  will  be  considered.  Radial 
elongation  of  the  arms  and  out-of-plane  deformations  are  neglected.  The  equations  of 
motion  for  the  uncontrolled  system  are  derived  by  finding  the  system  kinetic  and  potential 
energies,  using  assumed  modes  to  discretize  the  energies,  truncating  to  second-order  terms, 
forming  the  Lagrangian,  and  writing  Lagrange's  equations  for  the  system.  Figure  4-2 
shows  the  Draper/RPL  configuration  with  appendages  numbered,  deflections  and  rotations 
defined,  and  a  body-fixed  reference  frame. 
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Kinetic  Energy.  The  total  kinetic  energy  will  be  written  as  the  sum  of  the  hub, 
appendage  and  tip  mass  energies. 

The  kinetic  energy  of  the  rigid  hub  alone  is  that  of  a  solid  cylinder: 


Th  =  ^  ih 


(4-1) 


where  0  is  the  rigid-body  rotation  angle  of  the  body-fixed  reference  frame  with  respect  to 
an  inertial  reference  frame,  the  hub  moment  of  inerda  is  I),  =  m^  /  2,  mj,  is  the  hub  mass 
and  R  is  the  hub  radius. 

To  determine  the  kinetic  energy  of  appendage  1,  a  position  vector  is  written  as 


n  =  xb,  -h  vb2 


(4-2) 


The  velocities  of  the  elements  of  appendage  1  are  obtained  by  taking  the  first  time  derivative 
of  the  position  vector  with  respect  to  an  inertial  reference  frame: 


Ii  =  V  b2  +  (9x  £i) 


=  -  V  9  b]  -t-  (v  -t-  X  0)  b2 


(4-3) 


The  X  denotes  the  outer  or  cross  product.  Finally,  the  kinetic  energy  of  appendage  1  is  the 
inner  or  dot  product  of  the  velocity  vector  with  itself  integrated  with  respect  to  the  mass 
over  the  length  of  the  appendage.  Denoting  the  appendage’s  moment  of  inertia  as 


X  =  R  +  L 


X  dm 


(4-4) 
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the  kinetic  energy  of  appendage  1  is 
X  =  R  +  L 

Ti  =  -j  j  (li  •  Ii  )  dm 
x  =  R 


X  =  R  +  L 

M 

•• 

< 

1  (v^-*-v^0^-t-2v0x)  dm 

II  0^  > 

J 

L  x=  R  J 

(4-5) 


The  kinetic  energy  of  the  tip  mass  on  appendage  1  is  found  in  a  similar  manner.  The 
position  vector  is 


r  ti  -  (R  +  L)  bi  +  (v  I  X  =  R  +  L  )  ^2 


(4-6) 


As  a  reminder  that  the  deflection  v  must  be  evaluated  at  x  =  R  +  L,  let  Vjj  =  vl^  =  r  +  l  and 
write  the  velocity  vector  as 


I  tl  “  ''ll  ^2  ^  i  il^ 

=  —  Vjj  0  hi  +  (vti  +  0  (R  +  L)]  b2 


(4-7) 


The  velocity  is  then 

x/2 

Vti  =  III  •  Eti 

=  Vii  +  2  Vii  0  (R  +  L)  +  0  ^  (R  +  L)^  +  vj  0  ^  (4-8) 

To  include  the  small  effects  of  the  rotary  inertia  of  the  tip  mass,  the  angular  velocity  of  the 
tip  mass  is  that  of  the  body-fixed  reference  frame  plus  a  component  due  to  the  deflection  of 
the  end  of  the  appendage  (where  the  tip  mass  is  located): 
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(4-9) 


0.1  =  0  + 


0  +  v^i 


The  kinetic  energy  of  the  tip  mass  is  the  sum  of  the  translational  and  rotational  energies: 


(4-10) 


T,1  =  i 


The  kinetic  energy  of  appendage  3  and  its  tip  mass  are  derived  similarly: 


'  r  y  =  R  +  L  ] 

T3=  j  (u^  +  u^0^  +  2u0y)  dm  +  I3  0^  •  (4-11) 


Tg=  T  tmtvj  +  iiejl 


(4-12) 


where 


y  =  R  +  L 


h  = 


J  y^  dm 


Vj3  =  U(3  +  2  Uj3  0  (R  +  L)  +  0  ^  (R  +  L)  +  U(3  0  ^ 


(4-1 3a) 


(4- 13b) 


6ui3 


0,-1  =  0  + 


=  0  +  Ui3 


Ut3  =  Ul  y=R+L 


(4- 13c) 


(4- 13d) 
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Consider  only  the  case  of  anti-symmetric  deformations  of  the  arms  where  U4  =  -U3 
and  V2  =  -V2^  as  shown  in  Figure  4-3. 


Figure  4-3.  Anti-symmetric  deformation 

Thus,  the  kinetic  energies  of  the  tirst  (third)  appendage  and  tip  mass  are  equal  in  magnitude 
to  those  of  the  second  (fourth)  appendages  and  tip  mass  -  i.e.,  T2  =  Tj,  T(2  =  Tji,  T4  =  T3, 
Tt4  =  Tj3.  The  total  system  kinetic  energy  (T)  can  be  written  as  the  sum  of  the  individual 
energies: 

T  =  Th  +  2  ( Tj  Tti  +  T3  Tt3  )  (4-14) 
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To  expand  Eq  (4-14),  carry  out  the  integrations  over  dm  in  the  arm  thickness  (x)  and  height 
(K)  dimensions,  assume  constant  mass  density  along  the  arms  (pjj),  and  noting  that  I3  =  Ii, 
define 

lo  =  ^  Ih  +  2I1  -i-  2mi(R-t-L)^  (4-15) 


The  total  kinetic  energy  is  then  written  as 


T 


(Iq  +  2Ij)0^  -I-  p(jXK 


R  +  L 

J  (v^  -J-  e  ^  -I-  2  V  0  x)  dx 

R 


R  +  L 

+  J  (u^  +  u^  0^  2  u  0  y)  dy 

R 


I[  [  2  0  (vjj  -I-  U[3)  -t-  (vj])  -I-  (0^3)  ] 


-I-  mj  [  +  U{3  2  0  (R  -t-  L)  (v^j  +  Uj3)  +  0  ^  (vj^  +  uj)  ]  (4—16) 

A  discretized  system  model  is  formed  by  assuming  that  the  elastic  arm  deformations 
(relative  to  a  body-fixed  undeformed  state)  can  be  represented  as  a  linear  combination  of 
admissible  (shape)  functions  Pj  (p)  for  the  clamped-free  appendage  (21:699;  22:18): 

Pi(p)  =  1  -  cos(iTCp/L)  0.5(-l)i+l  (i7tp/L)2  (4-17) 

where  n  is  the  familiar  mathematical  constant,  p  is  the  radial  distance  along  the  arm 
(measured  from  the  end  clamped  at  the  hub)  and  L  is  the  arm  length.  This  assumed  mode 
or  shape  function  must  be  an  admissible  function  for  the  model  of  a  continuous  beam  with 
a  lumped  mass  at  the  end.  The  (shape)  functions  pj  (p)  are  admissible  since  they  satisfy 
the  geometric  boundary  conditions  of  the  appendage:  p  (p=0)  =  0  since  u  (t,0)  =  v  (t,0)  = 
0  and  p'  (p=0)  =  0  since  u'  (t,0)  =  v'  (t,0)  =  0.  [While  not  required,  the  Pj  (p)  are  also 
comparison  functions  for  a  clamped-free  appendage  (beam  alone),  since  they  also  satisfy 
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the  natural  or  physical  boundary  conditions:  p,"  (p=L)  =  0  since  u"  (t,L)  =  v"  (t,L)  =  0 
and  p'"  (p=L)  =  0  since  u"'  (t,L)  =  v'"  (t,L)  =  0.  However,  when  the  tip  masses  are 
included  in  the  model,  the  Pj  (p)  do  not  satisfy  the  natural  boundary  conditions  but  are  still 
admissible  functions.] 

The  transverse  body-fixed  deformations  of  the  arms  are  modeled  as 

Nsf 

U  (t,  p)  =  21  Hi  (p) .  0  <  p  <  L  (4-18a) 

i=  1 

Nsf 

v(t,  p)  =  2^  Vj  (t)  Pj  (p) ,  0<p<L  (4-18b) 

i=  1 

where  Njf  is  the  number  of  admissible  functions  considered.  This  assumed-modes  method 
yields  a  generalized  parameter  model  with  U,;  (t)  and  Vj  (t)  the  generalized  displacement 
coordinates.  Eqs  (4-18)  can  be  differentiated  with  respect  to  time  and/or  direction  and 
substituted  into  the  expression  for  the  total  system  kinetic  energy,  Eq  (4-16),  to  transform 
the  expression  to  generalized  coordinates. 

Potential  Energy.  The  potential  energy  of  the  system  is  the  sum  of  the  gravitational 
and  strain  energies.  Since  the  structure  is  assumed  to  be  in  orbit  (or,  if  on  earth,  supported 
by  an  air-bearing  table),  the  gravitational  potential  will  be  neglected  and  the  total  potential 
energy  P  taken  as  the  strain  or  elastic  potential  energy  alone,  summed  over  all  four 
appendages: 


Ey  is  the  modulus  of  elasticity  (Young's  modulus),  is  the  area  moment  of  inertia  (Ig  = 
kx3/12),  and  u"  and  v"  are  the  second  partial  derivatives  of  u(t,p)  and  v(t,p)  with  respect 
to  p.  The  restriction  of  anti-symmetric  deformations  also  remains. 

The  assumed-modes  method  [Eqs  (4-17)  to  (4-18)]  again  leads  to  the  generalized 
parameter  model  with  Uj  (t)  and  Vj  (t)  the  generalized  displacement  coordinates  and  ji,  (p) 
the  admissible  functions.  Substituting  into  generalized  parameters,  differentiating  and 
arranging  in  matrix  notation  yields 

P  =  Eyia  H,''(p)  |ij'(p)  dp  {U} 


Pi  (p)  Mj  (p) 


dp  {V} 


(4-20) 


System  Matrices.  Now  that  the  system  kinetic  and  potential  energies  are  in  hand,  the 
differential  equations  of  motion  are  obtained  from  Lagrange's  equation  (for  no 
nonconservative  forces) 


(4-21) 


where  i  =  1,2, ...,  #  of  generalized  coordinates.  The  Lagrangian  is  formed  as  Lg  =  T  -  P 
where  terms  of  order  three  or  higher  are  neglected  and  the  energies  are  written  as 

T  =  y  w"^  M  w  (4-22a) 


5wi, 

_  _ 

dt 

5  Wj 

p  =  1  ^ 


(4-22b) 
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M  and  K  are  the  mass  and  stiffness  matrices,  respectively  [as  shown  in  Eq  (3-1)],  and  the 
vector  of  generalized  coordinates  is 


yv  =  [  0  Ui  ,  U2  ,  -,  :  Vi  ,  V2  ,  -,  (4-23) 

Assuming  anti- symmetric  appendage  motion  has  reduced  the  number  of  generalized 
coordinates  by  2xNsf.  The  order  of  the  system  of  Eq  (3-1)  is  thus  n  =  (2xNsf)  +  1. 

The  elements  of  the  mass  and  stiffness  matrices  are,  for  i  =  1,2, Ngf  and  j  =  1,2, 

...,  Ngf 


M(l,l)  =  2(Io  +  It)  (4-24a) 

M(l,i-f-l)  =  2  [  m2  (R  +  L)  Pi  (L)  -h  It  |ij  (L) 

+  pd  X  K  f  (p  -h  R)  Pi  (p)  dp  ]  (4-24b) 

•'0 

M  (i+I,  j+1)  =  2  [  m2  Pi  (L)  pj  (L)  +  It  Pi  (L)  Pj  (L) 

-L 

+  Pd K  Pi  (p)  Pj  (p)  dp  ]  (4-24c) 

Jo 

M  (1,  i+Nsf-t-1)  =  M  (i+Nsf+1,  1)  =  M  (i+1,  1)  =  M  (1,  i+l)  (4-24d) 
M  (i+Nsfi-l,  j-HNsffl)  =  M  (i-hl,  j+1)  (4-24e) 

M(i+l,j+Nsffl)  =  M  (i+Ns{+l,j+l)  =  0  (4-24f) 
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K  (i+1,  j+1) 


r  »9  If 

2  Ey  la  Hi  (p)  pj  (p)  dp 


(4-25a) 


K(i+Nsffl,j+Nsf+1)  =  K(i+l.j+l)  (4-25b) 

K  (1,1)  =  K  (1,  i+1)  =  K  (1,  i+Nsf+1)  =  K  (i+1,  1) 

=  K(i+Nsf+1,  1)  =  K(i+l,j+Nsffl) 

=  K(i+Nsffl,j+1)  =  0  (4-25c) 


Model  Order 

After  the  first  (rigid-body)  mode,  conservation  of  angular  momentum  results  in  the 
vibrational  modes  occuring  as  pairs  of  'opposition'  and  'unison'  modes.  The  first  five 
normal  modes  are  shown  in  Figure  4-4.  The  anti-symmetric  'opposition'  modes  are  simple 
cantilever  beam  modes  where  adjacent  arms  move  in  opposition  such  that  equal  and 
opposite  torques  result  in  no  hub  (rigid-body)  rotation.  The  anti- symmetric  'unison' 
modes  are  characterized  by  all  four  arms  moving  in  unison;  consequently,  the  hub  rotates  to 
conserve  system  angular  momentum. 

The  frequencies  of  pairs  of  opposition  and  unison  modes  are  closely  spaced  and 
aligned  with  the  corresponding  frequency  of  simple  cantilever  beam  vibrations.  This 
alignment  and  decrease  in  spacing  with  increasing  mode  number  is  shown  in  Table  4-1  for 
the  first  nine  normal  modes  of  the  system  (one  rigid-body  mode  and  four  pairs  of 
opposition/unison  modes  aligned  with  the  first  four  cantilever  beam  modes). 
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Figure  4-4.  First  five  normal  modes  (21 JOO) 


i  able  4-1 .  Frequencies  for  first  nine  normal  modes 
Draper/RPL  Structure  Cantilever  Beam 


1 

Mode 

Frequency 

Mode 

Frequency 

# 

(Rad/sec) 

(Hz) 

# 

fRad/sec) 

(Hz) 

1 

0 

0 

• 

2 

4.40 

0.70 

3 

7.92 

1.26 

1 

9.42 

1.50 

4 

51.46 

8.19 

• 

5 

52.78 

8.40 

2 

59.04 

9.40 

6 

157.58 

25.08 

7 

158.34 

25.20 

3 

165.33 

26.31 

8 

313.09 

49.83 

# 

9 

314.35 

50.03 

4 

324.02 

51.57 

A  typical  design  approach  would  begin  with  a  very  large-dimensioned  structuial 
model,  perhaps  on  the  order  of  100  modes  or  more,  even  when  control  of  only  a  small 
number  of  lower  modes  is  desired.  Calculating  and  storing  eigenvalues  and  eigenvectors 
(frequencies  and  mode  shapes)  for  such  a  huge  number  of  modes  helps  ensure  the  accuracy 
of  the  model  at  the  lower  modes  by  Rayleigh-Ritz  convergence  to  the  true  mode  shapes  and 
frequencies  from  above.  The  model  for  control  design  is  then  a  very  accurate  but  truncated 
version  of  the  much  larger  structural  model.  However,  the  desire  in  this  research  to 
integrate  structural  and  control  design  precluded  such  an  approach,  as  the  basic  structural 
parameters  are  to  be  design  variables  and  iterated  each  loop  along  with  the  control  system 
variables.  Hence,  a  large-order  structural  model  was  not  calculated  and  truncated  to  form  a 
design  model.  The  structural  model  will  only  include  those  modes  to  be  controlled. 
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This  research  considers  control  system  design  by  specifying  damping  ratios  and 
frequencies  of  the  first  five  normal  modes,  i.e.,  the  rigid-body  mode  and  the  first  two  pairs 
of  vibrational  modes.  Thus,  N^f  =  2  and  w  is  a  five-dimensional  vector: 

w  =  [0,  U,,  U2.  V,,  V2]T  (4-26) 

Structural  system  design  is  accomplished  by  varying  the  thickness  and  height  of  the 
arms  and  hence  structural  mass.  The  arms  are  kept  identical  to  each  other  as  their 
dimension^  are  changed.  The  proportional  damping  constants  in  Eq  (3-4)  are  taken  to  be  a 
=  0,  p  =  10‘5,  so  C  =  10'^  K  and  damping  is  proportional  only  to  stiffness.  With  the 
values  for  the  system  parameters  given  in  Table  4-2,  the  explicit  expressions  for  the 
elements  of  the  M  and  K  matrices  [Eqs  (4-24)  and  (4-25)]  can  be  evaluated. 


Table  4-2.  Draper/RPL  configuration  parameters 


Hub  radius  (R) 

0.3048  m 

(1  ft) 

Rotary  inertia  of  hub  (1^) 

10.8465  Kg-m2 

(8  slug-ft2) 

Mass  density  of  arms  (p^) 

2,690  Kg/m^ 

(5.22  slug/ft3) 

Modulus  of  elasticity  (Ey) 

7.584xl0>0pa 

(1.584xl09lb/ft2) 

Arm  thickness  (x) 

0.003175  m 

(0.0104166  ft) 

Arm  height  (k) 

0.1524  m 

(0.5  ft) 

Arm  length  (L) 

1.2192  m 

(4.0  ft) 

Tip  mass  (mj) 

2.29038  Kg 

(0.156941  slug) 

Rotary  inertia  of  tip  mass  (Ij) 

0.00244  Kg-m2 

(0.0018  slug-ft2) 

Actuators  and  Sensors 

Torque  actuators  are  located  on  the  central  hub  and  at  some  position  py  on  each  arm. 
For  a  torque  uj,  at  the  hub,  a  torque  uj  at  position  Py|  on  arms  1  and  2,  and  a  torque  U2  at 
position  Pu2  on  arms  3  and  4,  the  right-hand  side  of  Eq  (3-1)  is 


Bu  =  Q  2ji’(pui) 


(4-27) 


2  U'  (Pu2)  J  ^  “2 


where 


iA'(p)  =  —  [U(p)].  JA(p)  =  tP]  (p).  (4-28) 

dp 


Angular  position  and  angular  velocity  sensors  are  located  on  the  hub,  and  colocated 
position  and  velocity  sensors  are  placed  at  stations  psi,  Ps2>  Ps3»  Ps4  3  and  4  and 

stations  ps5,  pj^,  ps7,  pgg  on  arms  1  and  2.  The  sensor  influence  matrix  H  is 


0  u’^(Psi) 


H  =  0  u'^(ps4) 


0  q'^  /(Pss) 


0  q'^  (Ps8) 


(4-29) 
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For  the  above  configuration,  then,  n  =  5  (Njf  =  2),  a  =  3  (actuators),  s  =  9  (sensors), 
B  is  a  5x3  (nxa)  matrix,  and  H  is  a  9x5  (sxn)  matrix.  The  gains  Gp  and  Gy  are  both  3x9 
(axs)  matrices,  and  the  output  vectors  are 

iLp  =  [  e,  u  (t,psi),  ...,  u  (t,ps4),  V  (t,ps5), ....  V  (t,ps8)  ]  T  (4-30a) 

Xv  =  Xp  (4-30b) 

The  locations  Puj  and  pu2  of  the  torque  actuators  and  pg],  P52.  Ps8  of  ^he  sensors 
will  be  allowed  to  vary  during  the  design  iterations  to  satisfy  design  objectives  of  placing 
closed-loop  frequencies  and  damping  ratios  (o)  and  Q  while  achieving  desired  reductions  in 
mass  and/or  control  effort.  Given  the  form  Eq  (3-3)  of  the  control  m.  the  control  effort  o  is 
defined  as  the  Frobenius  norm  of  the  10x10  (2nx2n)  matrix  G^u  given  by 


since 


'UU 


H'^Gj'GpH  h’^gJGvH 
Gv  Gp  H  h’’'  Gy  Gy  H 


(4-31) 


U^  U  =  JcT  X 


(4-32) 
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Note  that  if  a  solution  is  achieved  at  an  =  1,  in  multiobjective  optimization  terms  the 
solution  is  dominated  and  not  Pareto.  Also  notice  that  the  second  term  on  the  left-hand  side 
of  Eq  (5-3)  (term  in  { ))  may  be  interpreted  as  a  'search  direction'  multiplied  by  a  'step  size' 
(ttH).  This  formulation  has  favorable  convergence  characteristics  since  iterations  can  be 
started  with  close  estimates  of  the  solution  vector  [z(aH)]  at  each  continuation  step  as  is 
swept  from  zero  to  one.  The  'step  size'  (increment  in  ajj)  can  also  be  changed  adaptively 
in  each  iteration,  an  attractive  feature  for  highly  non-linear  objective  functions. 

For  a  given  continuation  step,  Eq  (5-3)  is  usually  underdetermined.  A  unique 
correction  vector  Az  is  obtained  by  minimizing  Az^  Az  subject  to  the  truncated  Taylor 
series  expansion  of  Eq  (5-3), 


Kz  +  Ag)  =  i(i)  + 


5z 


Ai  =  0 


(5-4) 


where  Wjj  is  an  arbitrary  weighting  matrix  and  i(z)  is  defined  as 

Ji  llCttH)]  =  ttH  Yi*  (i)  +  (1  -  ttH)  Yi  (Sstart)  *  Yi  i  =  1,  2 . p  (5-5) 

Conforming  to  the  implicit  local  linearity  assumptions,  the  weighted  minimum  norm 
correction  vector  is 


Az  =  -Wh' 


6z 


5z 


Wh' 


§z 


-1 


I(z) 


(5-6) 


where  the  pxd  locally  evaluated  Jacobian  is 


5z 

5z 

(5-7) 


Starting  with  a  neighboring  solution  at  the  previous  step  and  using  the  above 
correction  vector  Az,  z(aH)  is  refined  recursively  by 
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z(OtH)new  — (^•H^old 


(5-8) 


until  local  convergence  is  achieved  [Eq  (5-3)  is  satisfied]  for  each  ajj.  This  is  an  'inner 
loop'  on  Az.  Final  convergence  is  obtained  (an  'outer  loop')  by  incrementing  after  each 
local  convergence  until  an  =  1.  Table  5-1  summarizes  the  procedure. 

This  MCH  algorithm  has  proven  effective  for  eigenvalue  placement  problems  and,  by 
its  inherent  reliance  upon  neighboring  solutions,  is  also  attractive  for  systematically 
conducting  trade  studies  or  generating  trade-off  surfaces  (18-21,  23).  Junkins  (24)  also 
showed  that  this  minimum  norm  correction  technique  is  theoretically  equivalent  to  the 
gradient  projection  formulation  for  constrained  optimization  problems. 


Start 

Iteration  at 
—start  =  2  (0) 

2(aHi) 


2(aHq-l) 


Table  5-1.  MCH  procedure 

Iterate  Minimum 
Norm  Correction 

I[z(aHi)]  =  0 
I[z(aH2)]  =  0 


•a[z(l)]  =  0 


with  0  <  a^i  <  an2  <  .  .  ■  <  aHq=  1 


To  Obtain 
z(aHi) 

z(aH2) 


z(aHq=  1) 
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The  structural  design  objective  is  taken  to  be  the  mass  of  the  arms,  7i,  defined  as  the 
product  of  their  thickness,  height,  length  and  mass  density.  Control  design  objectives  are 
the  closed-loop  frequencies  co,  and  damping  ratios  Ci  of  the  first  five  modes  (one  rigid- body 
rotation  and  four  vibration  modes)  as  well  as  the  control  effort  o.  The  composite  objective 
vector  is 


=  [  (Di  ,  Cl  ,  0)2  ,  C2  >  -  .  “5  •  Cs  •  ^  ^  (5-9) 

The  design  variables  are  the  arm  thickness,  height,  actuator  and  sensor  locations,  and 
elements  of  the  feedback  gain  matrices.  The  design  vector  z  is 

z  =  [  "t,  K,  Pul,  Pu2’  Psl’  Ps2’  Ps8’  Gp(l*2),  "■»  Gp(l>9),  Gp(2,l), 

-,  Gp(3,9),  Gv(l,l),  Gy(l,2),  •••,  Gv(l,9),  Gv(2,l),  G^O,9)  ]  T  (5-10) 


The  arm  length,  L,  and  tip  mass,  mj,  have  also  been  included  as  design  variables  in 
previous  studies  (21, 22).  They  were  not  included  as  such  in  this  research  for  two  reasons: 
1)  they  are  more  likely  to  be  specified  due  to  higher-level  mission  requirements,  such  as 
keeping  some  equipment  away  from  the  hub,  while  the  arm  thickness  and  height  are  more 
likely  to  be  available  as  design  options  to  support  the  tip  masses,  and  2)  excluding  arm 
length  and  tip  mass  reduced  the  dimensionality  of  the  problem  and  the  consequent 
computational  burden  without  detracting  from  the  utility  and  validity  of  the  design  method. 

Results  of  studies  with  various  combinations  of  the  elements  of  2C  will  be  presented  in 
the  next  chapter. 

E)erivatives 

As  a  matter  of  implementation,  the  full  12x66  (pxd)  matrix  [Sj/Sj]  was  used  in  all 
design  iterations,  even  though  all  12  7j*  were  not  specified.  The  appropriate  Jj  were  set  to 
zero  in  Eq  (5-6)  to  calculate  correction  vectors  Az.  TTiis  was  acceptable  since  dominated 
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solutions  (tty  =  1)  were  sought  vice  Pareto  solutions,  which  seek  zero  gradients  to 
maximize  a^.  While  more  derivatives  were  calculated  than  necessary,  all  were  obtained 
analytically  (closed-form)  vice  numerically.  If  numerical  differentiation  had  been  required, 
the  additional  computational  burden  would  have  made  it  worthwhile  to  incorporate  the  logic 
to  manipulate  the  matrices  according  to  the  specific  objective  vector  under  consideration. 

To  detail  the  derivatives,  first  expand  the  right-hand  side  of  Eq  (5-7)  where  the  Yj  are 
as  given  in  Eq  (5-9): 


5Yi 

Syi 

ill 

5  Z| 

8  Z2 

8zd 

5Y2 

5Y2 

5Y2 

5  zj 

8  Z2 

8zd 

gyp 

.. 

^Yp 

8  zj 

8z2 

8zd 

(5-11) 


Since  tc  =  pj  t  k  L,  the  required  derivatives  are  simply 


5Yii  571  ,  Syu  57t  ,  /c  ^ 

-  =  -  =  Pd  K  L ,  -  =  -  =  pd  t  L  (5-12a) 

5zi  6x  5z2  6  k 


=  0  ,  i  =  2,  3,  •••,  d  =  66 


(5-12b) 


For  the  rates  of  change  of  the  closed-loop  frequencies  and  damping  ratios,  recall  Eq 
(3-13)  and  let  z^  be  any  element  of  the  design  vector  z; 
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6X  _  5a  ^  .  5cd 
5  5  Z|^  5  Z|^ 


(5-1 3a) 


(I)  a 


^Szic  j  iSz}^ 


To  obtain  the  5X/6zjj,  differentiation  of  Eq  (3-12a)  leads  to,  for  i  =  1, 2, 2n, 

S  ^  ^cl 

s  ~  R 

0  O 


(5-13b) 


(5-14) 


Given  the  form  of  A^i  in  Eq  (3-10), 


-1  (  5  Kci 
-M  ^  - 

V  6zj^ 


-1  f  ^  Ql 


(5-15a) 


-M  ‘  -  M  * 


(5- 15b) 


[Kgi  and  Qi  were  defined  in  Eq  (3-6)].  M  and  K  are  functions  of  x  and  k  only  (zj  and  Z2). 
Recalling  the  symmetry  of  those  matrices  given  in  Eqs  (4-24)  and  (4-25),  the  expressions 
for  the  non- zero  derivatives  are,  for  i  =  1,  2, ....  Njf  and  j  =1,2, Njf, 
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6M(1,1) 

1  Pd  K  Ii 

(5- 16a) 

5z] 

• 

5M(I,i-H) 

5  zj 

2  pd  K  f  (p  -t-  R)  (p)  dp 

Jo 

(5-1 6b) 

• 

5M(i-H,j-H) 

5zi 

2  pd  K  [  Pi  (p)  Pj  (p)  dp 

JO 

(5-16c) 

• 

5K(i-H,j-(-l) 

S  Zi 

1  Hj  (p)  p,'  (p)  dp 

(5-16d) 

5M(1,1) 

T  Pd  ^  II 

(5-1 7a) 

• 

5z2 

5M(l,i-(-l) 

5  Z2 

2pdX  f  (p-t-R)  Pi(p)  dp 

Jo 

(5- 17b) 

• 

6M(i-t-l,j+l) 

5z2 

2  Pd  “C  Pi  (P)  Pj  (P)  dp 

(5-17c) 

• 

5K(i-t-l,j-»-l) 

6  Z'l 

-L 

^EyX^I^  Pi  (p)  Mj  (p)  dp 

(5-1 7d) 

The  rest  of  the  components  of  and  Cd  are  not  functions  of  t  or  k.  Recalling  the 
definition  of  the  design  vector  z  [Eq  (5-10)], 


B  =  B  (Pui,  Pu2)  ” 

B  (Z3,  Z4) 

(5-1 8a) 

H  =  H  (Psi,  Ps2»  •••’ 

’  PsS)  ~  ^  •••'  ^12) 

(5- 18b) 

Gp  —  Gp  (Zi3'  ^14’  •• 

M  Z39) 

(5- 18c) 

Gy  —  Gy  (Z4Q,  Z4],  .. 

(5-18d) 
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Given  the  forms  of  B  and  H  [Eqs  (4-27)  to  (4-29)],  the  only  non-zero  elements  of  [5B/5z)j] 
and  [5H/5zjj]  are  in  the  positions  corresponding  to  the  particular  z^: 


5B 

—  =  2  m  (zk) ,  k  =  3, 4 

Szk 


(5-19a) 


5H  . 

7 —  =  Hi  (zk) . 

OZk 


k  =  5,  6 . 12 


(5- 19b) 


-  in 

Hi(p)  =  T7 


(5-l9c) 


Hi(p)  = 


i  TC  ^ 


/  ,xi+l  I  iTtp 

(-1)  +  COS  1  “L'y 


(5-19d) 


Note  i  =  1,  2  (i  =  1,  2, Njf;  Njf  =  2)  and  7t  is  the  mathematical  constant,  not  the  arm 
mass,  in  Eqs  (5-19). 

Even  more  simply,  the  elements  of  [6Gp/5zk]  and  [5Gv/5zk]  are  all  zeroes  except  for 
a  one  in  the  position  corresponding  to  the  particular  z^. 

Finally,  recalling  the  form  of  the  control  effort  matrix  Gyu  [Eq  (4-31)]  and  denoting 
the  magnitude  or  Frobenius  norm  of  G^^  as  IIGyu'!, 


\)  = 


1g 


uu ' 


r  2n  =  10  2n  =  10 
i=l  j=l 


(5-20a) 


c.  10  10 

5zk 


1=1  i=  1 


^  ^uu  ^^0^ 
5zk 


(5-20b) 


5-8 


Given  the  functional  dependencies  in  Eqs  (5-18)  and  the  expressions  for  [SH/Szj^]  in  Eqs 
(5- 19b)  and  (5- 19c),  the  elements  of  [SCuu/Sz^]  are 
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5G 


uu 


5zk 


'5Gv^ 

H 

<  8  zjj  y 

h't 

'8Gv' 

T 

GpH 

'5Gv" 

T 

Gy  +  G^ 

'5  Gy" 

}h 

. 

^  5  Zjj  y 

L  1 

^  8  zj.  ^ 

''  8  Z|j  y 

1  J 

k  =  40,41, ...,  66 


(5-2  Ic) 


Design  Procedure 

Now  that  all  necessary  equations  have  been  detailed,  the  following  are  the  steps  to  go 
through  one  'inner  loop'  of  the  design  procedure  (Table  5-1),  starting  with  ajj  =  ctni  = 
0.1  . 


1.  Initialization  -  set  constants,  read  initial  variable  values  from  input  files,  etc. 

2.  Calculate  system  matrices  -  M  [Eqs  (4-24)],  K  [Eqs  (4-25)],  B  [Eqs  (4-27), 
(4-28),  (5-19C)],  C  [Eq  (3-4)],  H  [Eqs  (4-29),  (4-17)],  [Eqs  (3-10),  (3-6)]. 

2..  Calculate  right  and  left  eigenvalues  and  eigenvectors  of  A^i  [Eq  (3- 1 1)];  sort  in 
order  of  increasing  frequency  and  normalize  eigenvectors  [Eq  (3-12)]. 

4.  Calculate  objective  vector  y  [Eq  (5-9)]  -  frequencies  and  damping  ratios  given 
by  Eq  (3-13),  control  effort  o  by  Eqs  (5-20a)  and  (4-31),  and  mass  7t  =  pj  t  k  L. 

1.  Calculate  gradients. 

g.  Partial  derivatives  of  mass  and  stiffness  matrices  [Eqs  (5-16)  and  (5-17)]. 

l2.  Partial  derivatives  of  control  and  sensor  influence  matrices  [Eqs  (5-19)]. 

£.  Partials  of  A^i  [Eqs  (5- 1 5),  (3-6)]. 
d.  Partials  of  closed-loop  eigenvalues  [Eq  (5-14)]. 

£.  Partials  of  closed-loop  frequencies  and  damping  ratios  [Eqs  (5-13)]. 
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f.  Partial  derivatives  of  mass  [Eqs  (5-12)]. 

g.  Partials  of  control  effort  [Eqs  (5-20),  (5-21),  (4-31)]. 

6.  Form  gradient  matrix  [5j/5^  [Eqs  (5-7),  (5-11)]  from  results  of  steps  (5-e)  thru 
(5-g)  above. 

2-  Calculate  correction  to  design  vector  [Eqs  (5-6),  (5-5)]. 

S..  Check  that  corrections  Az  will  not  violate  any  constraints  on  elements  of  z  (for 
example:  thickness,  height,  sensor  and  actuator  locations).  If  any  will  be  violated,  augment 
i  with  the  amount(s)  of  the  violation(s),  augment  [5J/5z]  with  a  -1  in  the  appropriate 
location(s)  [row(s)  and  column(s)],  and  re-calculate  correction  vector  Az  as  in  step  (7) 
above. 


9.  Calculate  magnitude  of  correction  vector  ( I  Az  1;  square  root  of  the  sum  of  the 
squares  of  the  elements). 

10.  If  I  Az  1  <  Z(qi  (convergence  limit),  this  inner  loop  has  converged  [Eq  (5-3)  is 
satisfied  at  this  value  of  a^].  Increment  an  and  begin  new  inner  loop  at  step  (7)  above. 

11.  If  inner  loop  has  not  converged,  calculate  new  design  vector  z  [Eq  (5-8)]  and 
return  to  step  (2)  above  for  another  iteration. 

Once  the  inner  loop  has  converged  [step  (10)]  at  an  =  1.  the  homotopy  family  of 
equations  [Eq  (5-3)]  has  become  the  original  system  of  Eq  (5-1),  and  the  procedure  has 
converged  to  the  desired  objective  vector  'x*. 
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Two  cases  were  considered,  most  readily  characterized  by  the  first  mode  eigenvalues. 
Case  1  is  high  frequency,  low  damping.  Case  2  is  low  frequency,  high  damping.  The 
desired  tOj*  and  were  not  specified  for  all  five  modes.  The  eigenvalue  portions  of  the 


objective  vectors  for  the  two  cases  are 

Case  1  :  7i*  =  [  cOj*  =  3  rad/sec,  Ci*  =  0.03,  ^2*  ”  0.03,  ^3*  =  0.03, 

C^*=0.01,  C5*=0.01]T  (6-1) 

Case  2:  Y2*  =  t  “1*  =  0-3,  0J2*  =  4.5.  CO3*  =  8.3,  Cl*  =  0.7, 

C2*  =  0.03,  Cs*  =  0.03,  C4*  =  0.01 ,  C5*  =  0.01  1  T  (6-2) 

Note  that  this  is  not  pole  placement  but  rather  partial  specification  of  the  eigenvalues  for  the 
first  fi  ,’e  normal  modes  of  vibration. 

All  elements  of  the  gain  matrices  Gp  and  Gy  were  initially  set  to  zero  except  for 
Gp(l,l),  which  was  set  to  O.OOi  to  prevent  the  rigid-body  eigenvalue  being  identically 
zero.  Initial  actuator  locations  were  Pyj  =  Pu2  =  L/2.  Initial  sensor  locations  were  p^j  = 
Ps5  =  L/4,  Ps2  =  Ps6  =  L/2,  Ps3  =  P57  =  0.7L,  Ps4  =  Ps8  =  0.9L.  Structural  parameters 
were  as  given  in  Table  4-2.  For  Case  2,  the  arm  length  and  tip  mass  were  adjusted  to  1.12 
m  (3.67  ft)  and  2.90  Kg  (0.19871  slug),  respectively,  to  place  the  second-  and  third-mode 
frequencies  (21). 

The  weighting  matrix  Wj^  in  Eq  (5-6)  was  nominally  set  to  the  identity  matrix. 
Typical  step  size  and  convergence  limits  were  Aa^  =  0.1,  Az  =  0.0001.  However,  the 
step  size,  weights  and  limits  on  the  sensor  and  actuator  locations  were  varied  to  aid 
convergence,  particularly  for  large  mass  and  control  effort  reductions. 
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Four  types  of  designs  were  accomplished:  1)  (partial)  eigenvalue  placement  [yi*  and 
72*  as  given  in  Eqs  (6-1)  and  (6-2)],  2)  eigenvalue  placement  with  mass  reduction,  3) 
eigenvalue  placement  with  control  effort  reduction,  and  4)  eigenvalue  placement  with  mass 
and  control  effon  reduction.  The  same  (partial)  eigenvalue  placement  was  used  for  all 
designs  within  each  case.  The  designs  are  identified  by  case  number  (first  digit)  and 
amount  of  mass  and/or  control  effort  reduction.  For  example,  design  1-Base  is  Case  1, 
baseline  (no  mass  or  control  effort  reduction),  design  1-MlO  is  Case  1,  10%  mass 
reduction,  design  2-CE25  is  Case  2,  25%  control  effort  reduction,  and  design  2-MCE50  is 
Case  2,  50%  mass  and  control  effort  reduction. 

Final  values  for  elements  of  the  design  vector  z  are  given  in  Appendices  A  (T,  K,  p^], 
Pu2’  Psb  Psl>  — ’  Ps8  designs)  and  B  (elements  of  Gp  and  for  designs  1-Base, 
1-M50,  1-CE50,  1-MCE50  and  2-MCE50).  Appendix  C  details  the  weighting  matrices 
Wpj  for  all  designs. 

Eigenvalue  Placement 

While  closed-loop  frequencies  and  damping  ratios  are  certainly  not  complete 
specifications  of  a  control  system,  they  are  familiar,  easily  understood  characteristics, 
particularly  for  top-level  concerns  of  structure/controller  interaction.  Solely  placing  the 
frequencies  and  damping  ratios  [as  given  in  Eqs  (6-1)  and  (6-2)]  was  not  the  purpose  of 
this  research,  but  is  common  to  all  integrated  structural/control  designs.  Designs  1-Base 
and  2-Base  were  accomplished  for  Cases  1  and  2,  respectively,  and  the  resulting  control 
efforts  u  were  calculated  [Eqs  (4-31)  and  (5-20a))  and  taken  as  the  baseline  values  [Case  1: 
12,5(K);  Case  2;  85].  Note  that  while  X)  does  not  have  a  physical  interpretation,  such  as 
torque,  it  is  still  a  valid  measure  of  the  amount  of  control  exerted  by  the  system  and,  hence, 
relative  comparisons  of  v  are  meaningful.  Actual  control  torques  require  specification  of 
the  state  vector. 
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Figures  6-1  and  6-2  are  typical  trajectories  of  the  eigenvalues  (Cases  1  and  2, 
respectively)  during  the  optimization  process  as  was  swept  from  =  0  (z  =  Zstart)  to 
Uh  =  1  (z  =  z*).  These  trajectories  are  representative  of  those  for  all  of  the  designs. 


Mass  Reduction 

For  integrated  stmctural/control  design,  the  previous  (partial)  eigenvalue  placement  is 
accomplished  simultaneously  with  reduction  of  the  arm  mass,  n.  The  objective  vectors  yi* 
and  Y2*  were  augmented  with  the  desired  arm  mass,  Jt*.  Mass  reductions  of  10,  25  and 
50%  were  obtained  while  placing  the  closed-loop  eigenvalues  as  specified  in  Eqs  (6-1)  and 
(6-2).  Mass  reductions  beyond  50%  were  not  attempted  but  are  certainly  possible  up  to 
that  corresponding  to  minimum  limits  on  the  arm  dimensions,  albeit  at  the  expense  of 
higher  control  efforts.  Control  efforts  for  the  reduced-mass  designs  were  calculated  and 
are  tabulated  in  Table  6-1.  As  expected,  maintaining  the  closed-loop  characteristics  while 
reducing  arm  mass  generally  requires  more  control  effort,  although  design  1-M25  (Case  1, 
25%  mass  reduction)  is  a  fortunate  exception,  apparently  due  to  more  favorable  sensor 
placements  (Appendix  A). 


Control  Effort  Reduction 

As  a  precursor  to  reducing  both  mass  and  control  effort,  control  effort  was  reduced 
while  holding  the  arm  mass  fixed.  The  objective  vectors  y\*  and  72*  were  augmented  with 
the  desired  control  effort,  \)*,  and  control  effort  reductions  of  10,  25  and  50%  of  the 
appropriate  baselines  were  obtained  while  placing  the  closed-loop  eigenvalues  as  specified 
in  Eqs  (6-1)  and  (6-2).  The  reductions  were  accomplished  primarily  through  sensor 
placement  (Appendix  A).  Reductions  beyond  50%  were  not  attempted. 
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Figure  6- 1 .  Case  1  eigenvalue  trajectories 
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Figure  6-2.  Case  2  eigenvalue  trajectories 


6- 


Table  6-1.  Control  efforts  for  reduced-mass  designs 


Design 

%  Mass 

Control  Effort 

# 

Reduction 

1-Base 

12,500 

1-MlO 

10 

13,500 

1-M25 

25 

11,000 

1-M50 

50 

32,400 

2-Base 

85 

2-MlO 

10 

339 

2-M25 

25 

1,441 

2-M50 

50 

4,323 

Mass  and  Control  Effort  Reduction 

Since  the  reduced-mass  designs  essentially  assumed  infinite  available  control  effort,  a 
more  realistic  problem  is  to  augment  the  basic  Yj*  and  Y2*  with  the  desired  arm  mass  snd 
desired  control  effort.  For  Case  1 ,  reductions  of  10,  25  and  50%  in  both  mass  and  control 
effort  were  obtained  while  placing  the  eigenvalues  per  Eq  (6-1). 

For  Case  2,  the  arm  length  and  tip  mass  had  been  adjusted  to  place  the  second-  and 
third-mode  frequencies  (21)  and  the  first-mode  frequency  is  very  low  (0.3  rad/sec). 
Consequently,  the  baseline  control  effort  (required  to  satisfy  Y2*)  is  also  very  low,  so  low 
that  only  5%  reductions  in  mass  and  control  effort  could  simultaneously  be  achieved 
(design  2-MCE5).  Control  efforts  much  larger  than  the  baseline  were  required  to  satisfy 
eigenvalue  placement  with  simultaneous  mass  reduction,  as  shown  in  Table  6-1.  These 
control  efforts  were  subsequently  reduced  by  the  same  percentage  as  the  mass,  i.e.,  the 
control  efforts  corresponding  to  10,  25  and  50%  mass  reductions  were  reduced  by  10,  25 
and  50%,  respectively  (designs  2-MCElO,  2-MCE25,  2-MCE50). 
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Although  design  goals  of  mass  and/or  control  effort  were  met  while  achieving  desired 
closed-loop  frequencies  and  damping  ratios,  the  design  model  only  included  the  first  five 
normal  modes.  Model-order  tmncation,  while  a  necessary  evil  for  computational  purposes, 
leaves  the  designs  susceptible  to  spillover  from  the  higher-order  unmodeled  or  residual 
modes.  All  designs  achieved  previously  (Appendix  A)  were  evaluated  against  7-,  9-  and 
11-mode  evaluation  models.  Figures  6-3  through  6-7  are  plots  of  the  modes'  damping 
ratio  for  designs  1-Base,  1-M50,  1-CE50,  1-MCE50  and  2-MCE50  (Case  1:  baseline, 
50%  mass  reduction,  50%  control  effort  reduction,  50%  mass  and  control  effort  reduction, 
and  Case  2:  50%  mass  and  control  effort  reduction)  to  show  the  extreme  spillover 
sensitivity  of  the  designs.  The  modes  of  the  design  model  (1-5)  are  those  whose  damping 
ratios  were  specified  during  the  design  process.  The  'baseline'  damping  (solid  bar)  is:  1) 
the  specified  damping  ratio  per  Eq  (6-1)  for  modes  1-5,  and  2)  the  open-loop  damping 
ratios  for  the  residual  modes  (slighdy  stable  due  to  the  assumed  proportional  damping). 

All  the  designs  had  at  least  one  mode  go  unstable,  while  many  unstable  modes  was 
more  of  the  norm.  Instability  was  not  limited  to  the  residual  modes;  modes  in  the  design 
model  were  also  driven  unstable,  as  with  modes  2  and  5  in  Figure  6-6  and  mode  4  in 
Figure  6-7.  Even  if  the  modes  in  the  design  model  were  not  driven  unstable,  they  were 
often  subject  to  wild  variations  in  damping  ratio  as  soon  as  any  higher-order  modes  were 
included  in  the  evaluation  model  (modes  4  and  5  in  Figure  6-5,  mode  4  in  Figure  6-6). 
There  was  no  quantifiable  correlation  of  spillover  sensitivity  with  mass  and/or  control  effort 
reduction,  either  within  or  across  the  two  cases,  other  than  all  designs  were  very 
susceptible. 
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Figure  6-3.  Damping  ratio  spillover  sensitivity,  design  1-Base 
(Case  1 ,  baseline) 
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Figure  6-4.  Damping  ratio  spillover  sensitivity,  design  1-M50 
(Case  1 , 50%  mass  reduction) 
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Figure  6-5.  Damping  ratio  spillover  sensitivity,  design  l-CESO 
(Case  1, 50%  control  effort  reduction) 
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Figure  6-6.  Damping  ratio  spillover  sensitivity,  design  1-MCE50 
(Case  1 ,  50%  mass  and  control  effort  reduction) 
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Figure  6-7.  Damping  ratio  spillover  sensitivity,  design  2-MCE50 
(Case  2,  50%  mass  and  control  effort  reduction) 
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The  minimum  correction  homotopy  technique  described  in  Chapter  V  has 
successfully  been  used  for  integrated  structural/control  design  as  evidenced  by  the  designs 
achieved  (listed  in  Appendix  A).  However,  in  light  of  the  design  model  order  truncation 
and  lack  of  any  robustness  criteria  in  the  design  objective  vector,  the  extreme  spillover 
sensitivity  of  the  designs  is  not  surprising.  Although  spillover  does  not  have  to  be 
destabilizing,  it  will  certainly  affect  system  performance.  Thus,  in  order  to  provide  a 
complete  technique  for  integrated  structural/control  design,  spillover  control  is  also 
addressed.  Two  approaches  were  considered:  simple  constraints  and  modal  suppression. 

Constrained  Solutions 

The  first  approach  to  control  spillover  was  simply  to  increase  the  order  of  the  design 
model  from  5  to  1 1  modes.  No  attempt  was  made  to  place  the  eigenvalues  of  the  additional 
modes,  6-11;  they  were  simply  constrained  to  have  damping  ratios  >  0.(X)1.  In  effect,  this 
was  an  attempt  to  simply  meet  the  design  objectives  given  in  Chapter  VI  while  guaranteeing 
stability  through  1 1  modes;  i.e.,  push  off  the  onset  of  spillover  to  modes  12  and  beyond, 
hoping  that  these  much  higher-frequency  modes  have  less  effect  on  the  modes  of  the 
original  design  model  (1-5).  Or,  it  may  be  thought  of  as  a  'brute  force'  attempt  to  build  a 
deadband  into  the  controller  over  the  frequency  range  covered  by  modes  6-11  (~  52-840 
rad/sec). 

As  a  matter  of  implementation,  the  objective  vector  was  expanded  to  include  the  first 
1 1  modes 


^  = 


[  a>i  ,  ^ I  ,  0)2  ,  C2  >  ■"  >  1  >  Cl  1  >  ^  ^ 


(7-1) 
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[mass  (7t)  and  control  effon  (o)  were  included  as  appropriate].  Even  though  all  24  y*  were 
not  specified,  the  full  24x66  matrix  [5i/5z]  was  used  in  all  design  iterations.  The 
appropriate  Jj  were  set  to  zero  in  Eq  (5-6)  to  calculate  correction  vectors  Az.  Since  this 
approach  was  only  applied  to  the  Case  1  designs,  the  Jj  corresponding  to  0)2,  0)3, ....  Wn 
were  all  set  to  zero.  The  spillover  control  then  set  the  Jj  corresponding  to  i  =  6, 7, ..., 
11  equal  to  zero  if  that  ^j  was  >  0.001.  Finally,  design  iterations  were  started  with  the 
'converged'  solutions  from  the  Case  1  designs,  i.e.,  the  final  designs  were  'post- 
processed'  to  control  spillover. 

The  new  designs  (with  spillover  control)  were  then  evaluated  against  13-,  15-,  17-, 
19-  and  21 -mode  evaluarion  models.  Figures  7-1  and  7-2  are  plots  of  the  modes'  damping 
ratios  for  designs  1-Base  and  I-MCE50  (Case  1,  baseline  and  Case  1,  50%  mass  and 
control  effort  reduction)  with  the  spillover  control  included.  Figure  7-1  (design  1-Base) 
shows  the  design  stable  and  not  too  much  variation  in  damping  ratio  through  21  modes. 
However,  this  design  was  the  exception  rather  than  the  rule.  The  rest  of  the  spillover- 
control  designs  exhibited  characteristics  similar  to  the  non-spillover-control  designs:  all 
had  multiple  modes  go  unstable,  including  the  spillover-controlled  modes  (6-11),  and  the 
damping  ratios  were  often  subject  to  wild  variations  with  changing  evaluation  model  size. 
The  25  and  50%  mass  and  control  effort  reduction  designs  (1-MCE25, 1-MCE50)  also  had 
modes  in  the  design  model  driven  unstable  (modes  2  and  5  in  Figure  7-2). 

In  short,  the  simple  constraint  approach  to  spillover  control  did  delay  the  onset  of 
instability  to  modes  12  and  higher,  but  the  designs  still  exhibit  far  too  much  spillover 
sensitivity.  Model-order  truncation  effects  were  not  mitigated  by  this  technique  except  for 
one  fortunate  design  (1-Base),  and  that  very  exception  points  out  the  unreliability  of  this 
approach  to  control  spillover. 
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Figure  7-1.  Damping  ratio  spillover  sensitivity  with  simple  constraint  spillover  control, 

design  1-Base  (Case  1,  baseline) 
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Figure  7-2.  Damping  ratio  spillover  sensitivity  with  simple  constraint  spillover  control, 
•  design  1 -MCE50  (Case  1 , 50%  mass  and  control  effon  reduction) 
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A  modal  suppression  technique  (25, 26)  was  attempted  next.  While  the  technique  has 
also  been  used  to  design  decoupled  controllers  (27-29),  it  is  used  here  to  eliminate 
observation  spillover  and  recover  stability. 

Technique.  Recalling  the  system  of  Eq  (3-1),  let  and  ifcj,  i  =  1,  2, ...,  n  be  the 
eigenvalues  and  eigenvectors  of  the  associated  eigenvalue  problem 

\|/i2  M  (|ii  =  K  iUi  (7-2) 

Define  O  as  the  matrix  with  as  the  /th  column,  where  the  iljj  have  been  normalized  with 
respect  to  the  mass  matrix, 

<1)T  M  O  =  1  (7-3) 


i.e.,  O  is  the  system  modal  matrix.  Also  define  the  coordinate  transformation 

w  =  <1>  n  (7-4) 

where  n  is  an  n-vector  of  modal  coordinates. 

The  second-order  system  of  Eq  (3-1)  can  now  be  written  in  first-order  state-space 
form  as 


X  =  Ax  4^  Dm 


where 


(7-5) 


(7-6) 


7-4 


A  = 


0  I 

L  r-2^vj 


(7-7) 


D  = 


0 


(7-8) 


r is  a  diagonal  matrix  of  open-loop  modal  frequencies  squared,  and  is  a 

diagonal  damping  matrix.  The  output  y  and  control  u  are 


X  = 


H<D  0 
0  Hd) 


(7-9) 


ii  =  [  Go  :  Gy  ]  y 


To  eliminate  observation  spillover,  first  define  a  new  output  q  as 


a  =  r  y 


(7-10) 


(7-11) 


The  matrix  F  will  be  chosen  such  that  q  does  not  contain  any  information  regarding  modes 
to  be  suppressed,  say  modes  6-r.  In  other  words,  q  will  be  a  linear  combination  of  the  s- 
elements  of  the  output  y  such  that  only  modes  1-5  are  observed.  Since  there  are  nine 
sensors  (s=9),  if  five  modes  are  to  be  controlled,  a  maximum  of  four  modes  can  be 
suppressed  (i.e.,  r=9,  suppress  modes  6-9).  There  must  be  as  many  independent  'sources 
of  information'  (elements  of  the  new  output  q)  as  modes  to  be  controlled. 

The  columns  of  HC>  correspond  to  observations  of  modes: 


H  d>  =  a  =  [  w]  :  £,2  '•  ’ ' '  •  J 


(7-12) 
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To  suppress  observations  of  modes  6-9,  F  must  be  such  that 

r  [  Sg  :  H7  :  Eg  :  29  ]  =  F  E5_9  =  0  (7-13) 

Singular  value  decomposition  is  used  to  determine  the  null  space  of  and  F  solved  for 
by  transposing  Eq  (7-13) 

E^9  =  0  (7-14) 

The  new  control  becomes 

u  =  [  Gp  :  Gy  ]  q  =  [  Gp  ;  Gy  ]  F  y 

FE  0 

X  (7-15) 

0  FE_ 

Note  that  the  gain  matrices  Gp  and  G^  must  be  reduced  in  size  accordingly.  For 
suppression  of  modes  6-9,  Gp  and  Gy  are  3x5  matrices. 

Results.  The  'converged'  solutions  from  eigenvalue  placement,  50%  mass  reduction, 
50%  control  effort  reduction,  and  simultaneous  50%  mass  and  control  effort  reduction 
were  'post-processed'  to  eliminate  observation  spillover;  1 )  structural  sizes  and  sensor  and 
actuator  locations  were  held  fixed,  2)  F  was  calculated  to  suppress  observation  of  modes 
6-9,  and  3)  the  MCH  algorithm  was  re-run  with  the  elements  of  the  3x5  Gp  and  Gy  gain 
matrices  the  only  variables  adjusted  to  recover  the  closed-loop  eigenvalue  placement. 
Control  effort  was  not  specified  (i.e.,  the  desired  objective  vectors  were  as  shown  in  Eqs 
(6-1)  and  (6-2),  unaugmented]  and  the  5-mode  design  model  was  used. 

Tables  7-1  through  7-3  detail  th"  damping  ratios  obtained  with  up  to  a  19-mode 
evaluation  model  for  designs  1-Base,  1-MCE50  and  2-MCE50  (Case  1,  baseline;  Case  1, 
50%  mass  and  conn’ol  effort  reduction;  Case  2,  50%  mass  and  control  effort  reduction). 


However,  all  modal-suppression  designs  exhibit  the  same  favorable  characteristics.  The 
controlled  modes'  (1-5)  damping  varies  only  slightly;  there  appears  to  be  some 
'conservation  of  damping  ratio'  within  mode  pairs  (2&3,  4&5)  at  higher-order  evaluations. 
The  suppressed  modes  (6-9)  remain  at  their  open-loop  damping  ratios,  as  they  should; 
since  the  controller  can't  'see'  those  modes,  it  doesn't  affect  them.  The  assumed  slight 
amount  of  proportional  damping  keeps  modes  6-9  stable.  Not  only  are  the  controlled  and 
suppressed  modes  much  more  well-behaved,  but  the  residual  modes'  (10-19)  damping 
ratios  are  relatively  immune  to  changes  in  the  order  of  the  evaluation  model.  There  is  now 
essentially  one  decade  of  deadband  (modes  6-9,  ~  52-530  rad/sec)  above  the  controller 
bandwidth  (modes  1-5, 0-52  rad/sec). 

Two  designs  were  also  attempted  suppressing  only  modes  6-7  (note  Gp  and  Gy  were 
3x7  matrices).  However,  modes  8-10  remained  unstable,  whereas  when  modes  6-9  were 
suppressed,  the  onset  of  spillover  instability  was  delayed  until  mode  1  *t  the  earliest  (for 
those  two  specific  designs).  As  one  might  expect,  making  full  use  of  the  four  extra  sensors 
(suppressing  modes  6-9)  gave  more  robust  designs  than  only  putting  in  two  modes  of 
deadband  (modes  6-7). 

Finally,  Table  7-3  shows  the  design  stable  through  19  modes.  This  fortunate  (as 
opposed  to  planned)  result  points  out  the  necessity  to  always  check  through  all  modes  of 
interest  to  see  if  residual  modes  are  stable  or  unstable  -  spillover  does  not  necessarily 
always  cause  instability. 


Table  7-1.  Damping  ratios  with  modal  suppression,  design  1-Base 

(Case  1,  baseline) 


#  of  Modes  in  Evaluation  Model 


Mod?# 

2 

11 

ii 

li 

12 

1 

0.03 

0.0299 

0.03 

0.0301 

0.0301 

2 

0.03 

0.0299 

0.03 

0.03 

0.03 

3 

0.03 

0.0301 

0.03 

0.03 

0.03 

4 

0.01 

0.0106 

0.0105 

0.0104 

0.0103 

5 

0.01 

0.00973 

0.00978 

0.00983 

0.00985 

6 

0.000799 

0.000799 

0.000799 

0.000799 

0.000799 

7 

0.000803 

0.000803 

0.000803 

0.000803 

0.000803 

8 

0.0016 

0.0016 

0.0016 

0.0016 

0.0016 

9 

0.00161 

0.00161 

0.00161 

0.00161 

0.00161 

10 

0.0324 

0.0327 

0.0323 

0.0322 

11 

0.0357 

0.0346 

0.0349 

0.035 

12 

0.0166 

0.0165 

0.0164 

13 

-0.05 

-0.0499 

-0.0499 

14 

-0.00609 

-0.00608 

15 

-0.00788 

-0.00786 

16 

-0.00314 

17 

0.00874 

•  18 
19 
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12 

0.0301 

0.03 

0.03 

0.0102 

0.00986 

0.000799 

0.000803 

0.0016 

0.00161 

0.032 

0.0351 

0.0164 

-0.0496 

-0.0066 

-0.00751 

-0.00314 

0.00875 

0.0287 

0.0141 


Table  7-2.  Damping  ratios  with  modal  suppression,  design  1-MCE50 
(Case  1 ,  50%  mass  and  control  effort  reduction) 


#  of  Modes  in  Evaluation  Model 


Mode# 

2 

11 

ii 

ii 

11 

12 

1 

0.03 

0.029 

0.029 

0.0292 

0.0295 

0.0295 

2 

0.03 

0.0299 

0.03 

0.0298 

0.0298 

0.0297 

3 

0.03 

0.0297 

0.0296 

0.0296 

0.0296 

0.0296 

4 

0.01 

0.0113 

0.0114 

0.0113 

0.0111 

0.0111 

5 

0.01 

0.00998 

0.00998 

0.00998 

0.00998 

0.00998 

6 

0.000867 

0.000867 

0.000867 

0.000867 

0.000867 

0.000867 

7 

0.000869 

0.000869 

0.000869 

0.000869 

0.000869 

0.000869 

8 

0.00171 

0.00171 

0.00171 

0.00171 

0.00171 

0.00171 

9 

0.00172 

0.00172 

0.00172 

0.00172 

0.00172 

0.00172 

10 

0.0755 

0.0756 

0.0755 

0.0753 

0.0753 

11 

0.00877 

0.00876 

0.00876 

0.00876 

0.00876 

12 

0.00609 

0.006 

0.00591 

0.00591 

13 

0.00694 

0.00703 

0.00706 

0.00708 

14 

0.00443 

0.00444 

0.00444 

15 

-0.00784 

-0.00747 

-0.00753 

16 

-0.035 

-0.0348 

17 

0.00924 

0.00924 

• 

• 

Table  7-3. 

Damping  ratios  with  modal  suppression,  design  2-MCE50 

(Case  2, 50%  mass  and  control  effort  reduction) 

#  of  Modes  in  Evaluation  Model 

• 

Mode# 

1 

u 

iA 

li 

12 

12 

1 

0.7 

0.699 

0.7 

0.7 

0.7 

0.7 

2 

0.03 

0.0297 

0.0298 

0.0298 

0.0298 

0.0298 

A 

3 

0.03 

0.0306 

0.0305 

0.0304 

0.0304 

0.0305 

w 

4 

0.01 

0.00602 

0.00642 

0.00691 

0.00657 

0.00646 

5 

0.01 

0.0138 

0.0134 

0.013 

0.0133 

0.0134 

6 

0.000958 

0.000958 

0.000958 

0.000958 

0.000958 

0.000958 

• 

7 

0.00096 

0.00096 

0.00096 

0.00096 

0.00096 

0.00096 

8 

0.00189 

0.00189 

0.00189 

0.00189 

0.00189 

0.00189 

9 

0.00189 

0.00189 

0.00189 

0.00189 

0.00189 

0.00189 

• 

10 

0.00404 

0.00406 

0.00401 

0.00407 

0.0041 

11 

0.00363 

0.00355 

0.00361 

0.00354 

0.00349 

12 

0.00554 

0.00555 

0.00557 

0.00558 

13 

0.00563 

0.00562 

0  00561 

0.0056 

• 

14 

0.007 

0.00709 

0.00712 

15 

0.00723 

0.00705 

0.00697 

16 

0.0119 

0.0121 

17 

0.0116 

0.0115 

• 

18 

0.0149 

19 

0.0151 

• 

• 
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Vni.  Conclusions 


This  dissertation  has  presented  an  approach  for  integrated  structural/control  design 
and  applied  it  to  a  flexible  space  structure  (Draper/RPL  configuration)  and  its  active  control 
system  (direct  output  feedback).  The  primary  contribution  has  been  to  treat  the  problem  in 
an  integrated  or  simultaneous  fashion  vice  sequential  design  iterations  on  the  two  separate 
problems. 

A  minimum  correction  homotopy  (MCH)  technique  was  used  to  simultaneously 
optimize  a  vector  of  objective  functions,  including  closed-loop  damping  ratios  and 
frequencies,  structural  mass,  and  control  effort.  Hence,  designs  were  achieved  via 
multiobjective  optimization  vice  forming  a  scalar  objective  via  weighting  and  summing  the 
individual  objectives.  The  MCH  technique  has  heretofore  not  been  used  to  solve 
multiobjective  optimization  problems,  nor  has  structural  design  (reduce  mass)  been 
combined  with  control  system  design  at  such  a  detailed  level  (place  closed-loop  frequencies 
and  damping  ratios,  reduce  control  effort).  [Note  that  the  control  system  design  was  not 
pole  placement,  but  specification  of  all  damping  ratios  and  some  frequencies  of  interest,  as 
shown  in  Eqs  (6-1)  and  (6-2)]. 

Designs  with  up  to  50%  reductions  in  mass  and/or  control  effort  were  obtained  while 
achieving  desired  c'osed-loop  characteristics  for  the  first  five  normal  modes  of  vibration. 
However,  the  designs,  based  on  a  reduced-order  structural  model,  are  easily  driven 
unstable  by  spillover  from  higher-order  unmodeled  modes.  A  modal  suppression 
technique  completed  the  design  approach  by  eliminating  observation  spillover  from  modes 
6-9  and  providing  one  decade  of  deadband  above  the  controller  bandwidth.  The  resultant 
designs  are  much  less  sensitive  to  the  effects  of  model-order  truncation.  Additional 
deadband  could  easily  be  obtained  by  controlling  fewer  modes  or  simply  adding  sensors. 


IX.  Recommendations 


The  designs  obtained  in  this  research  are  dominated  vice  Pareto  solutions.  That  is,  all 
elements  of  the  objective  vector  x  were  simultaneously  achievable;  no  trade-off  amongst 
the  elements  of  the  objective  vector  (yj,  individual  objectives)  was  done  to  arrive  at  some 
subjective  'best'  solution.  Attempts  to  use  both  utility  function  and  goal  programming 
techniques  as  described  in  Chapter  II  to  find  Pareto  optimal  solutions  to  this  multiobjective 
optimization  problem  were  unsuccessful.  In  fact,  mass  and  control  effort  were  never 
included  as  objectives,  since  convergence  for  the  sub-problem  of  eigenvalue  placement  was 
never  achieved.  A  variety  of  weights  and  goals  were  tried,  not  only  with  closed-loop 
frequencies  and  damping  ratios  as  design  objectives  but  also  using  the  real  and  imaginary 
parts  of  the  closed-loop  eigenvalues  directly  as  the  objectives  in  an  attempt  to  at  least  reduce 
the  nonlinearities  of  the  problem,  but  all  to  no  avail. 

The  first,  natural  extension  of  this  research  would  be  to  use  the  minimum  correction 
homotopy  (MCH)  design  approach  to  find  Pareto  solutions.  The  reported  body  of  research 
on  multiobjective  optimization  and  Pareto  solutions,  albeit  not  for  integrated 
structural/control  design  but  rather  dealing  with  only  one  discipline,  suggests  utility 
function,  goal  programming,  and/or  game  theory  techniques  (2-11,  13-17).  Since  all 
reported  techniques  involve  some  form  of  weighting  and  summing  individual  design 
objectives,  direct  application  of  the  MCH  algorithm  may  be  difficult.  However,  given  the 
excellent  convergence  characteristics  of  the  MCH  algorithm  to  dominated  solutions,  it 
seems  worthwhile  to  attempt  to  use  it  to  find  Pareto  solutions. 

Whether  dominated  or  Pareto  solutions  are  sought,  a  valid  question  is;  what  are  the 
best  objective  functions  to  use  for  integrated  structural/control  design?  The  design 
objectives  used  in  this  research  (closed-loop  frequencies  and  damping  ratios,  structural 
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mass,  control  effort),  although  certainly  valid  and  important  system  characteristics,  can 
certainly  be  improved  upon.  Additional  objectives  could  be  eigenvectors,  robustness, 
structural  stresses,  displacements,  etc.  Any  design  improvements  would  then  have  to  be 
weighed  against  the  corresponding  computational  burdens  and  difficulties. 

Any  or  all  of  the  assumptions  and  parameters  inherent  in  this  proof-of-technique 
could  be  investigated.  The  design  approach  certainly  needs  to  be  examined  with  respect  to 
much  higher-dimensional  problems,  allthough  no  stumbling  blocks  other  than 
computational  loading  are  immediately  foreseen.  Using  the  technique  with  different 
structural  models  and/or  control  laws  should  not  only  verify  the  technique's  utility  and 
validity  but  also  yield  great  varieties  of  design  options  and  trade-offs. 

Finally,  the  modal  suppression  technique  could  be  incorporated  into  the  main  design 
iteration  loop  (vice  being  a  'post-processing'  algorithm  to  recover  stability)  to  see  if  such 
inclusion  achieves  required  stability  with  lower  control  gains,  control  effort  or  sensitivity  to 
uncertain  parameters. 
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A  -  1 


2-MCE25  0.01080  0.3618  0.945  0.950  0.418  0.466  0.491  0.626  0.411  0.433  0.458  0.748 

2-MCE50  0.01153  0.2259  0.980  0.980  0.371  0.385  0.374  0.509  0.361  0.378  0.379  0.527 
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Appendix  B:  Final  Design  Gains 
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Appendix  C:  Weiehting  Matrices 


All  weighting  matrices  Wjj  [in  Eq  (5-6)]  were  diagonal  matrices.  For  all  but  two 
designs,  the  matrices  were  identity  matrices. 

For  design  1-MCE50,  the  first  12  diagonal  elements  were  1,  1,  10,  10,  10,  10,  100, 
100,  10,  10,  100,  100;  the  rest  were  ones  -  i.  e.,  the  elements  of  corresponding  to  the 
actuator  locations  and  two  inboard  sensor  locations  on  each  arm  (p^j,  Py2.  Psi»  Ps2>  Ps5> 
Ps5)  were  10,  and  the  weights  on  the  two  outboard  sensors  on  each  arm  (P53,  Ps4,  pj-;, 
Psg)  were  100. 

For  design  2-CE50,  the  first  12  diagonal  elements  were  1,  1,  0.1,  C.i,  0.1,  0.1,  0.1, 
0.1,  0.1,  0.1,  0.1,  0.1;  the  rest  were  ones  -  i.e.,  the  elements  of  Wj3  corresponding  to  the 
sensor  and  actuator  locations  were  all  0.1. 
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